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Abstract 

r*"*- ■ A semi-microscopic model for nucleon pairing in nuclei is presented starting from the ab intio BCS 

m 

gap equation with Argonne Vxs force and the self-consistent Energy Density Functional Method 
fT} ' basis characterized with the bare nucleon mass. The BCS theory is formulated in terms of the 

o ' 

model space So with the effective pairing interaction calculated from the first principles in the 
subsidiary space S' . This effective interaction is supplemented with a small phenomenological 



^ \ addendum containing one phenomenological parameter universal for all medium and heavy atomic 

?H ■ 

nuclei. We consider the latter as a phenomenological way to take into account both the many-body 
corrections to the BCS theory and the effective mass effects. For protons, the Coulomb interaction 
is introduced directly. Calculations made for several isotopic and isotonic chains of semi-magic 
nuclei confirm the validity of the model. The role of the self-consistent basis is stressed. 
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I. INTRODUCTION 



In the last few years, some progress has been made in the microscopic theory of pairing 
in nuclei by the Milan group [l-3] and by Duguet et al. In the first paper of 

the Milan series, the Bardeen-Cooper-Shrieffer (BCS) gap equation for neutrons with the 
Argonne V14 potential was solved for the nucleus 120 Sn which, being located in the middle 
of the tin isotopic chain, is a traditional benchmark for the nuclear pairing problem. The 
Saxon- Woods Shell-Model basis with the bare neutron mass m* = m was used, and the 
discretization method in a spherical box was applied to take into account the continuum 
states restricted by the limiting energy -E max = 600 MeV. Rather optimistic results were 
obtained for the gap value, Abcs — 2.2 MeV. Although it is bigger of the experimental 



.eft the hope to achieve a good 
3j the basis was enlarged, i.e. 



one, A exp — 1.3 MeV, the difference is not so dramatic and 
agreement by developing corrections to the scheme. In Refs. 
-Emax = 800 MeV, and, what is more important, the effective mass m* 7^ m was introduced 
into the gap equation. The new basis was calculated within the Skyrme-Hartree-Fock (SHF) 
method with the Sly4 force p], that makes the effective mass m*(r) coordinate dependent, 
remarkably different from the bare one m. E.g., in nuclear matter the Sly4 effective mass is 
equal to m* = 0.7m. As it is well known, in the so called weak coupling limit of the BCS 
theory, the gap is exponentially dependent, i.e. A oc exp(l/g), on the inverse dimensionless 
pairing strength g = m*V e sk-F/ir 2 , where V e g is the effective pairing interaction. Therefore, 
a strong suppression of the gap takes place in the case of m* < m. The value of Abcs = 0.7 
MeV was obtained in Ref. [2J and Abcs = 1-04 MeV, in Ref. In both cases, the too 
small value of the gap was explained by invoking various many-body corrections to the BCS 
approximation. The main correction is due to the exchange of low-lying surface vibrations 
("phonons"), contributing to the gap around 0.7 MeV ffl, so that the sum turns out to be 
A = 1.4 MeV very close to the experimental value. In 3j the contribution of the induced 
interaction caused by exchange of the high-lying in-volume excitations was added either, 
the sum again is equal to A ~ 1.4 MeV. Thus, the calculations of Refs. 

aa 

showed that 

the effects of m* 7^ m and of many-body corrections to the BCS theory are necessary to 
explain the difference of (Abcs — A exp ). In addition, they are of different sign and partially 
compensate each other. Unfortunately, both effects contain large uncertainties. Indeed, the 
calculations 

aa 

were carried out using the Skyrme force parameters, that are fixed only 
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in the vicinity of the Fermi surface, whereas the BCS gap equation involves the function 
m*(k) for k 3> k?- The same is true for the induced interaction. In this case, the situation is 
even worse since the spin channel in the high-frequency response function is important, but 
the corresponding combinations of Skyrme force parameters are not known sufficiently well 
even at the Fermi surface. More detailed analysis of these problems can be found in Ref. 
Q. On the other hand, a similar compensation between self-energy and vertex corrections 
is found in the ca.cu.ation of pairing in nuclear mattergB, supported with a Monte Carlo 
calculation of the gap[10[. 

A "blow-up" was produced in 2008 by Duguet and Losinsky {4] who solved the ab initio 
BCS gap equation for a lot of nuclei on the same footing. It should be noticed that the main 
difficulty of the direct method of solving the nuclear pairing problem comes from the rather 
slow convergence of the sums over intermediate states A in the gap equation, because of the 
short-range of the free A A- force. Evidently, this is the reason why the authors of Refs. 2, ^] 
concentrated only on a single nucleus, i.e. 120 Sn. To avoid the slow convergence, the authors 
of Refs. 2, 5 1 used the "low-k" force Mow-k fl which is m fact very soft. It is defined 
in such a way that it describes correctly the AA-scattering phase shifts at momenta k<A, 
where A is a parameter corresponding to the limiting energy E\ [m ~ 300 MeV. Moreover, 
the force Vi ow _k rapidly vanishes for k>A, so that in the gap equation one can restrict 
the energy range to _E max ~300 MeV . In addition a separable version of this force was 
constructed that made it possible to calculate neutron and proton pairing gaps for a lot of 
nuclei. Usually the low-k force is found starting from some realistic AA-potential V with 
the help of the Renormalization Group method, and the result does not practically depend 
on the particular choice of V In addition, in Ref. [4| V[ ow _k was found starting from the 
Argonne potential vi§, that is different only a little bit from Argonne V14, used in Ref. sj]. 
Finally, in Ref. the same SLy4 self-consistent basis was used as in Ref. Thus, the 
inputs of the two calculations look very similar, but the results turned out to be strongly 
different. In fact, in Ref. jj] the value Abcs — 1-6 MeV was obtained for the same nucleus 
-So which is already b, gg er than the experimental one by * 0.3 MeV. In Refs. fl □ 
we analyzed the reasons of these contradictions. This point was discussed also in Ref. |5j. 
It turned out that, in fact, these two calculations differ in the way they take into account 
the effective mass. Namely, in Refs. {2, 3] the effective mass m* = m^ hy4 was used for all 
momenta up to k max = y/2m*E max ~ 6 fin -1 . On the other hand, in Refs. {4, 5] this mass 
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is used only for k < A ~ 3 fm _1 , and the prescription m* = m is, in fact, imposed for 
bigger momenta. It implies that the gap A depends not only on the value of the effective 
mass at the Fermi surface, as it follows from the above exponential formula for the gap, 
but also on the behavior of the function m*(k) in a wide momentum range. However, this 
quantity is not known sufficiently well [7|], which makes rather uncertain the predictions of 
such calculations. 

To avoid such a drawback, we suggest a semi-microscopic model for nuclear pairing con- 
taining a single phenomenological parameter. It starts from the ab initio BCS gap equation 
with the Argonne force v 18 treated with the two-step method. The complete Hilbert space S 
of the problem is split into the model subspace So of low-energy states and the complemen- 
tary one S'. The gap equation is solved in the model space with the effective interaction V e g 
which is found in the complementary subspace. This ab-initio term of V e g is supplemented 
by a small one-parameter addendum that should hopefully embody all corrections to the 
simplest BCS scheme with m* = m. Preliminary results of this model are reported in Ref. 

Q. 



II. OUTLINE OF THE FORMALISM 



We start from the general form of the many-body theory equation for the pairing gap A 

A T =U T G T G S T A T , (1) 

where r = (n, p) is the isotopic index, U T is the AA-interaction block irreducible in the two- 
particle r-channel, and G T (Gy) is the one-particle Green function without (with) pairing. 
A symbolic multiplication, as usual, denotes the integration over energy and intermediate 
coordinates and summation over spin variables as well. We have used above the term "BCS 
theory" meaning that, first, the block U of irreducible interaction diagrams is replaced with 
the free AA-potential V in Eq. (pQ), and, second, the simple quasi-particle Green functions 
for G and G s are used, i.e. those without phonon corrections and so on. In this case, Eq. ([1]) 
turns greatly simplified and can be reduced to the form usual in the Bogolyubov method, 

A r = -V T x T , (2) 
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where 

de 



is the anomalous density matrix which can be expressed explicitly in terms of the Bogolyubov 
functions u and v, 

x r (n,r 2 ) = ^<(ri)i)[(r 2 ). (4) 

i 

Summation in Eq. (Hj) scans the complete set of Bogolyubov functions with eigen-energies 
Ei > 0. 

As mentioned in the Introduction, we use a two-step renormalization method of solving 
the gap equation in nuclei to overcome the slow convergence problem. We split the complete 
Hilbert space of the pairing problem S to the model subspace So, including the single-particle 
states with energies less than a fixed value of E , and the subsidiary one, S'. The gap 
equation is solved in the model space: 

A r = V: s G T G s T A T \ So , (5) 

with the effective pairing interaction V^ s instead of the block U T in the original gap equation 
(JTJ). It obeys the Bethe-Goldstone type equation in the subsidiary space, 

v; s = u T + u T G T G T v; s \ S '. (6) 

In this equation, the pairing effects can be neglected provided the model space is sufficiently 
large. That is why we replaced the Green function G S T for the superfluid system with its 
counterpart G T for the normal system. In the BCS approximation, the block W in Eq. (jHj) 
should be replaced by V T . To solve this BCS version of Eq. (jSJ) in non- homogeneous systems, 
we have found a new form of the local approximation, the Local Potential Approximation 



(LPA). Originally, it was developed for semi- infinite nuclear matter 15 1, then for the slab 

n n □ n 

of nuclear matter (see review articles [16|, [17|) and, finally, for finite nuclei [7|, [12[. It turned 
out that, with a very high accuracy, at each value of the average cm. coordinate R = 
(ri + r2 + r3 + r4)/4, one can use in Eq. ([6]) the formulae of the infinite system embedded into 
the constant potential well U = U(R). This explains the term LPA, and also significantly 
simplifies the equation for V e g, in comparison with the initial equation for A. As a result, 
the subspace S' can be chosen as large as necessary to achieve the convergence. From 
the comparison of the direct solution of Eq. fl6]) in the slab with LPA, it was shown that 
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LPA has high accuracy, even in the surface region, for sufficiently large model space, Eq 
~20-=-30 MeV). For finite nuclei, including 120 Sn, the validity of LPA was also checked 
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121 ] . In this case, the boundary energy should be made larger up to E =40 MeV. In 



this paper, we use LPA with this choice of En for systematic calculations of the gap in 
spherical nuclei. For V, we use, just as in Ref. pf, the Argonne potential v 18 . To make the 
calculations more treatable, we use the separable representation [3] of the vi 8 potential. 
Even in this simplified version the calculations of the set of matrix elements of V e g for a 
single nucleus require about 30 hours cpu with 50 processors of the multi-processor system 
of the Kurchatov Institute. 

Let us notice that the use of the low-k force Vi ow _k could be also interpreted in terms 
of the two-step renormalization scheme of solving the gap equation ([T]), with E'q— 300 MeV 
and with free nucleon Green functions G in Eq. ([6]) (i.e. U(R) = 0). Then, (with U— >V) 
one obtains V e g— >V\ ow --k (see Ref. [19( where the usual renormalization scheme is used to 
find V[ ow _k instead of the Renormalization Group equation). Now, the comparison of the 
direct solution of the gap equation (CQ) (or (T5])) in Ref. |3|] with the Argonne iViV-potential 
V and of "renormalized" equation ([5]) with V e s = Vi ow -k shows that the difference appears 
because, in the subsidiary subspace S', the effective mass m*^m is used in the first case, 
and m* = m in the second one. Thus, the resulting gap depends not only on the value of 
the effective mass at the Fermi surface, but also on the behavior of the function m*(k) in a 
wide momentum range. This dependence was demonstrated explicitly in Refs. [3, Q]- The 
use of the SHF effective mass corresponding to the SLy4 force, or to any other version of the 
Skyrme force, could hardly be accepted. Indeed, these effective forces were introduced and 
fitted to describe systematically nuclear masses and radii. As a rule, the description of the 
single-particle spectrum nearby the Fermi surface with Skyrme forces is rather poor, and it 
is expected not to be less poor at those high momenta that are involved in the gap equation 
([!]). This point makes it tricky the problem of determine the pairing gap completely from 
first principles, because the many-body theory(CE]) contains, in addition to the "/c-mass" of 



the SHF method, the "E'-mass" (inverse Z-factor) [8|, |20j, l21[ , that is not sufficiently well 
known even in nuclear matter [7- 9] . The corrections to the BCS version of Eq. ([1]) include 
also the difference of the block U from the potential V, mainly due to the so-called induced 
interaction. The attempt in Ref. [3] to determine the latter from the SLy4 force together 
with the nuclear mean field looks questionable. Indeed, the SLy4 parameters were fitted to 
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the nuclear mass table data mainly related to the scalar Landau-Migdal (LM) amplitudes 
/, /'. As to the spin amplitudes g, g', they remain practically undetermined in the SHF 
method. But the contribution of the spin channel to the induced interaction is not smaller 
than that of the scalar one [3]. Parameters g,g' are well known from the calculations of 
nuclear magnetic moments within the Finite Fermi Systems (FFS) theory |22| but, as for 
the Skyrme parameters, only at the Fermi surface. However, the states distant from the 
Fermi surface are important to calculate the induced interaction. The induced interaction 
for such states has only been determined in nuclear matter within the microscopic Brueckner 
theory {9]. At last, let us imagine to get from phenomenology the functions m*(k), Z(k) and 
all the LM amplitudes far from the Fermi surface. Even in this case, the use of so many 
phenomenological ingredients devalues significantly the ab initio starting point, i.e. the free 
NN potential V in the pairing gap calculation. 

Instead, we suggest to introduce in the effective pairing interaction a small phenomeno- 
logical addendum which embodies, of course approximately, all the corrections to the BCS 
scheme discussed above. The simplest ansatz for it is as follows: 

V; ff ( ri , r 2 , r 3 , r 4 ) = V^{r u r 2 , r 3 , r 4 ) + ^Gq^S^ - v 2 )5{v 1 - r 3 )5(r 2 - r 4 ). (7) 

Here p(r) is the density of nucleons of the kind under consideration, and 7 1 " are dimensionless 
phenomenological parameters. To avoid any influence of the shell fluctuations in the value 
of p(0), the average central density p(0) is used in the denominator of the additional term. 
It is averaged over the interval of r<2 fm. The first, ab initio, term in the r.h.s. of Eq. (J7J) 
is the solution of the BCS version of Eq. (jfJJ) (with hi — > V) in the framework of the LPA 
method described above, with m*=m in the subspace S' . 

We will see below that a rather small value of the phenomenological parameter 7„ = 
7 P ~ 0.06 is sufficient to produce the necessary effect of suppressing theoretical gaps pre- 
dicted by the ab initio calculation. The smallness of the phenomenological addendum to 
the effective interaction itself is demonstrated in Fig. 1 where the localized "Fermi av- 
erage" effective interaction is drawn for 7 = and 7 = 0.06 values for two heavy nu- 
clei. In the mixed coordinate-momentum representation, this quantity is defined as follows: 
V ff(ki, k 2 , ri, r 2 ) -)> V^ S (R = ri)5(ri - r 2 )d(r 1 - r 3 )5(r 2 - r 4 ), where 

V e F ff (#) = J d 3 tV eS {h = k 2 = k F (R), R - t/2, R + t/2), (8) 
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FIG. 1: The Fermi average effective pairing interaction Vf s (R) for 120 Sn and 200 Pb nuclei, at 7 = 
(dashed curves) and 7 = 0.076 (solid curves) 

with kp(R) = i/2m(/i — U (R)), provided /i — U(R) > 0, and k F (R) = otherwise. Here 
fi and U(R) are the chemical potential and the potential well, respectively, of the kind of 
nucleons under consideration. A similar quantity was considered before in the nuclear slab 
to visualize the effective interaction properties 16|, |32|. At a glance, the difference between 
the interaction strengths for 7=0 and 7=0.07 is negligible, but it produces noticeable effects 
in the gap owing to the exponential dependence on the force of the gap discussed in the 
Introduction. 

In the case of proton pairing (r = p), we retain in Eq. ([6]) only the "strong" part of the 
pp potential V p and then add the "bare" Coulomb potential Vc to the BCS term of Eq. (J7J) 

= VfS + V C . (9) 

Such approximation is valid because the mixed strong- Coulomb term of Eq. fl6]) is short- 
range, just as the strong term itself, but it is proportional to the small parameter of the fine 
structure constant a = 1/137. For matrix elements (A1A2IVCIA3A4) of the bare Coulomb force 
this small parameter is partially compensated in the diagonal case, Ax — A2 — A3 — A4 — A, 
owing to its long-range character. Such matrix elements can be estimated as ~ e 2 /R, where 
R is the nuclear radius. For example, for N = 82 isotonic chain it is of the order of 15-20% 
of the main diagonal matrix elements. It results in a~ 30% suppression of the proton gap 
value. So strong contribution of the Coulomb interaction to the proton pairing was reported 
recently in Ref. jsj. The above arguments to neglect the renormalization of Coulomb 
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interaction are valid also inside the model space. In particular, the photon vertex T(q) does 
not change in the long-range (small q) limit, owing to the Ward identity. In other words, 
the large (long-range) Coulomb matrix elements are not changed by the strong interaction. 
Small (short-range) matrix elements do change, but they can be neglected. Thus, we can 
put 7 P = 7„ = 7 in Eq. (j7|) after adding the bare Coulomb interaction to the BCS term for 
protons. 

Then, the gap equation (jSJ) in the model space is solved in the A-representation, with the 



self-consistent 
method 



24 



jasis determined within the Generalized Energy Density Functional (GEDF) 



28], where m*=m is assumed, with the functional DF3 271, |28| . The latter 
is of principal importance for our approach: first, because it makes the results less model- 
dependent, all effects of m* ^ m in both model and subsidiary subspaces being attributed to 
the in-medium corrections beyond the pure BCS approximation, and second, because single- 
particle spectra of the GEDF method rule, in a better agreement with experiment 



than those of the popular versions of the SHF method [29( (see for comparison Ref. [30]). 
The quality of the single-particle spectrum nearby the Fermi surface is very important for 
obtaining the correct value of the gap calculated from Eq. (j2J). 

III. RESULTS 

We solved the equations of the previous Section using the self-consistent A-basis of the 
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. The discretization method for 



GEDF method with the DF3 functional of Refs. 
the continuum states was used in the spherical box of radius R = 16 fm with the grid 
step h = 0.05 fm. The model space So was extended up to the energy E = 40 MeV, the 
subsidiary one 5", up to E max = 1000 MeV. The numerical stability of the results was checked 
by increasing the parameters up to Eq = 60 MeV, E max = 1200 MeV and R = 24 fm, and 
we found for the gap value a numerical accuracy of 0.01 MeV. 

In this paper, we limit ourselves to semi-magic nuclei. There are several reasons for 
such a choice. The first one is of technical nature, namely we have only spherical code 
for the gap equation, and all or almost all semi-magic nuclei are spherical. The second 
one is that in nuclei with both non-magic subsystems there are often very "soft" low-lying 
2 + -states whose contributions to different quantities, in particular to the gap value, can be 
rather strong and non-regular. In such cases, it is difficult to expect that the simple ansatz 
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(J7J) will work with an universal parameter 7. For neutron pairing, we consider the lead, 
tin, and calcium chains. The formulae above correspond to so-called "developed pairing" 
approximation [ijjl, that amounts to impose the equality of the A + and A~ operators and to 
neglect the particle-number non-conservation effects. Therefore we limit ourselves to nuclei 
having, as a minimum, four particles (holes) above (below) the magic core. For this reason, 
only the isotope 44 Ca was considered in the calcium chain. 

In accordance with the recipe of Ref . Q , we represent the theoretical gap with the "Fermi 
average" combination 

A F = ^(2j+l)A AA /E^ +1 )' ( 10 ) 

A A 

where the summation is carried out over the states A in the interval of \e\—fj,\<3 MeV, . 
The "experimental" gap is determined by the usual 5-term mass difference, Eq. flA3j) of 
the Appendix. As it is argued in the Appendix, the relevance of the mass difference to the 
gap has an accuracy of ~ (0.1 4- 0.2) MeV. Therefore, it is reasonable to try to achieve the 
agreement of the gap within such limits. 

Let us begin from the lead chain. Results are presented in Table I. In the end of this 
table, the results are shown for the 44 Ca nucleus which is the lightest one among the scope 
"from calcium to lead" considered in this article. We try to find a value of 7 which will be 
universal for all the region. 

In Table I the strong effect of the small phenomenological addendum in Eq. ([7]) is shown. 
The ab initio BCS result (7 = 0) significantly overestimates the gap. Switching on this term 
with 7 = 0.06 4- 0.08 suppresses the gap by 30 - 40%, in agreement with the data. The 
rms deviation of the theoretical values of the gap from the data for 13 nuclei, presented in 



Table I, is W(5A) 2 ~0.045 MeV at 7 = 0.06 and 0.103 MeV at 7 = 0.08. It has a minimum 



(5A) 2 ~0.037 MeV at 7 = 0.064, but, according to the above estimate, it is not reasonable 
to push the accuracy too much. In any case, we may consider the parameter 7 ~ 0.06 as an 
optimal one for this set of nuclei. 

Let us consider now the tin chain. The results are presented in Table II. Again we see a 
strong suppression with 7 = 0.06 4- 0.08, but the agreement now is remarkably poorer than 



in the lead case. Now, the rms deviation is \ (<5A) 2 ~0.165 MeV at 7 = 0.06 and 0.169 at 



7 = 0.08, and the minimal value y (<5A) 2 ~0.158 MeV at 7 = 0.07 is also too large. 

Let us move now to protons. The effect of the Coulomb interaction to the proton gap is 
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TABLE I: Neutron gap Ap (MeV) in Pb isotopes and 44 Ca nucleus. 



nucleus 


7=0 


An 
F 

0.06 


0.08 


^cxp 


182 Pb 


1.79 


1.33 


1. 


.20 


1.30 


184 Pb 


1.79 


1.33 


1. 


.20 


1.34 


186p b 


1.78 


1.32 


1. 


.19 


1.30 


188p b 


1.76 


1.31 


1. 


.17 


1.25 


190p b 


1.73 


1.29 


1. 


.16 


1.24 


192p b 


1.68 


1.22 


1. 


.09 


1.21 


194p b 


1.62 


1.16 


1. 


.03 


1.13 


196p b 


1.53 


1.09 


0. 


.96 


1.01 


198p b 


1.43 


1.00 


0. 


,87 


0.94 


200 p b 


1.31 


0.90 


0. 


.80 


0.87 


202p b 


1.16 


0.79 


0. 


.69 


0.78 


204p b 


0.95 


0.64 


0. 


.56 


0.71 


44 Ca 


1.83 


1.50 


1.41 


1.54 



shown in Table III for the isotonic chain of nuclei with the magic neutron number N = 82. It 
is seen that, indeed, it is rather strong ~ 0.5 MeV, in accordance with jlj]. Again at 7=0.06 
the agreement is almost perfect for the most part of nuclei, and only for the two heaviest 
isotones the disagreement is of the order of 0.2 MeV. In this case, a possible explanation 
lies in that we are close to the phase transition to a deformed state (at A ~ 150). Owing 
to the contribution of these "bad" cases, the average difference between the theoretical and 
experimental gaps for the = 82 chain is rather high, J (8 A) 2 = 0.124 MeV. The average 
error for all 34 nuclei considered is equal to J (<5A) 2 ~0.086 MeV. As it follows from the 
analysis discussed in the Appendix, this value is within the accuracy of the experimental 
values of the gap extracted from the 5-term formula flA3j) . 
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TABLE II: Neutron gap A£ (MeV) in Sn isotopes. 



nucleus Ap A, 



7 =0 0.06 0.08 



106 Sn 


1 


35 


0.95 





83 


1 


20 


108 Sn 


1 


52 


1.13 


1 


01 


1 


23 


110 Sn 


1 


65 


1.26 


1 


14 


1 


30 


112 Sn 


1 


74 


1.34 


1 


23 


1 


29 


114 Sn 


1 


80 


1.40 


1 


28 


1 


14 


116 Sn 


1 


82 


1.43 


1 


31 


1 


10 


118 Sn 


1 


83 


1.44 


1 


32 


1 


25 


120 Sn 


1 


80 


1.42 


1 


31 


1 


32 


122 Sn 


1 


74 


1.38 


1 


28 


1 


30 


124 Sn 


1 


65 


1.30 


1 


21 


1 


25 


126 Sn 


1 


51 


1.19 


1 


10 


1 


20 


128 Sn 


1 


31 


1.02 





94 


1 


16 



IV. THE ROLE OF THE SINGLE-PARTICLE SPECTRUM 



We consider the poor agreement for the tin chain as a troubling point. Indeed, this chain 
is a traditional benchmark for the pairing problem in nuclei, and the most strong deviation 
takes place for the 116 Sn nucleus which is in the very center of the chain where the scheme 
used should work especially well. In searching the reasons for such a drawback, we paid 



attention to a new version of the DF3 functional, named DF3a 



30] , in which the spin-orbit 



and effective tensor terms of the original DF3 functional were changed to fit new data on 



spin-orbit splitting in magic nuclei 



311 ] . The matter is that the details of neutron paring 



in the chain under consideration depend essentially on the position of the "intruder" state 
l/in/2 which, in turn, is determined mainly by these components of the functional. We 
repeated the calculations for this functional DF3a and found that the agreement became 



much better. The optimal value 7 = 0.052 with y (<5A) 2 ~0.056 MeV is now a little less than 
for the lead chain, but for 7 = 0.06 chosen above the agreement is also reasonably good, 



(5A) 2 ~0.084 MeV. In Fig. 2 is displayed the comparison of the theoretical predictions for 
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TABLE III: Proton gap A p (MeV) for the isotone gap N = 82. 



nucleus Ap A exp 

v c p ff =v c ° ff V c p ff =V c ff + Vc 







7=0 


0.06 




136 Xe 


1.65 


1.19 


0.87 


0.75 


138 Ba 


1.80 


1.33 


0.98 


0.87 
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FIG. 2: Neutron gap in Sn isotopes 

the pairing gap in the tin chain (both versions of the functional under consideration) with 
the experimental data. 

We see that, indeed, the new calculation is now in nice agreement with the data. To 
find the reason of so strong difference of the results for these versions of the essentially same 
functional, we examined the corresponding single-particle spectra along the chain comparing 
them with existing experimental data. Dealing with an even A Sn isotope, there is a dilemma 
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FIG. 3: The /in/2 ^vel position accounted from the ground state in Sn isotopes 

how to consider a state \i) under consideration, either the hole state (i.e. the excitation of 
the A_1 Sn nucleus) or the particle one (the excitation of the A+1 Sn nucleus). We use a 
simple recipe: the state \i) is considered as a hole state if the inequality vf > 0.5 takes place 
and as a particle state otherwise. Note that in the case of vf ~ 0.5 the difference between 
the particle and hole energies is, as a rule, quite small. In general, both the functionals 
reproduce the experimental low-lying levels sufficiently well. In particular, the spin of the 
ground state of odd isotopes is always reproduced correctly for both calculations. But there 
is a noticeable difference for the l/in/2 state in the left part of the chain, till 122 Sn, see Fig. 
3. 

In this region, its position is systematically lower than the experimental one for the DF3 
functional and systematically higher, for the DF3a functional. In the vicinity of the 116 Sn 
nucleus, this difference becomes dramatic and the too low position of the l/in/2 state in the 
DF3 case leads to an incorrect enhancement of the gap value. 

In Fig. 4, the single-particle spectrum of the 114 Sn nucleus calculated for the two versions 
of the functional is compared with the experimental one. Each theoretical level is supple- 
mented with the (2ji + l)uiVi factor that determines mainly the contribution of the i-level 
to the gap equation. The same for 116 Sn nucleus is displayed in Fig. 5. We see that for all 
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FIG. 5: 116 Sn spectrum 

other states \i) these factors are rather close for the two versions of the functional, but for 
the lhn/2 state the difference is rather strong. As to heavy tin isotopes for which this level 
becomes the ground state for both functionals, the difference between their predictions for 
the gap value are quite close, see Fig. 2. 

This analysis shows the great sensitivity of the gap value to the single-particle spectrum 
nearby the Fermi level, especially to the position of levels with high j -value. Therefore it 
is interesting to examine which effect is to be expected in the other cases of going from the 
initial DF3 functional to this new version. It is displayed in Fig. 6 for the lead isotopes. 
We see that in this case the overall agreement for the new version of the functional becomes 
worse. Evidently, again the position of high j-levels is different, in favor of the DF3 functional 
in this case. 

A bad situation for the DF3a functional rises for the N = 82 isotone chain, see Fig. 7. 
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FIG. 6: Neutron gap in Pb isotopes. 
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FIG. 7: Proton gap for N = 82 isotones 

The analysis shows that again the lh n / 2 level is guilty, now for protons. For the DF3a 
functional it is again higher than the experimental position, and for 144 Sm and 146 Gd nuclei 
much higher. As the result, it does not practically contribute to the gap equation resulting 
in a too small gap value. 
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V. CONCLUSION 



We suggest a simple semi-microscopic model for the nuclear pairing starting from the 
ab initio BCS gap equation. The self-consistent GEDF Method basis, characterized by the 
bare nucleon mass, is employed in the calculation. The gap equation is recast in the model 
space So, replacing the bare interaction with the effective pairing interaction determined in 
the complementary subspace 5". 

The Argonne vxs potential was adopted along with the LPA method. A small phenomeno- 
logical term is added to this effective interaction that contains one parameter which should 
embody approximately the effective mass and various corrections to the pure BCS theory. 



Calculations are carried out with the DF3 functional [27|, |28| for semi-magic lead and tin 
isotopic chains and the N = 82 isotonic chain as well. The Coulomb interaction is explic- 
itly included in the proton gap equation. We find that the model reproduces rather well 
the experimental values of the neutron and proton gaps in semi-magic nuclei. The overall 
agreement (a/ (<5A) 2 ~0.13 MeV) is better than that obtained in Ref. [4j], where the authors 
did not introduce free parameters explicitly but they made it implicitly by using a specific 
/c-dependence in the effective mass. 

We examine also the role of the single-particle spectrum in the gap equation. For this 



aim we use the new modification DF3a [3fJ of the functional |27|, |28| that changes spin-orbit 
and effective tensor terms. The use of this functional gives an agreement better for the tin 
chain and worse for the lead chain and even more for the N = 82 chain. The accuracy of 
the predictions depends strongly on the quality of reproducing the positions of high j-levels 
in the self-consistent basis used. We are thinking, e.g., to the l/in/2 neutron level in the tin 
isotopes and proton level for N = 82 isotones. 

The ansatz of Eq. flTJ) exhibits an obvious drawback. The phenomenological GEDF 



pairing interaction of Ref. 



28] contains the surface term (<x(dp/dr) 2 ) that plays an essential 



role for the description of the odd-even effect (staggering) in nuclear radii. It originates 
mainly from the exchange by surface phonons which was explicitly taken into account in 

a a. 

The addition, to such a term in Eq. (JJJ) is associated the introduction of a new 
parameter, and at the first stage we prefer to avoid that. A more consistent scheme should, 
evidently, include the explicit consideration of the low-lying phonons, as e.g. in [2|, but 
taking into account the so-called tadpole diagrams 33]. In this case, the phenomenological 
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constant 7 should, of course, change. 
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Appendix A: Accuracy of extracting experimental gap values from the mass differ- 
ences 

In this Appendix, we discuss the accuracy of determination of the "experimental" gap, 
A exp , from the mass data. Usually, this quantity is found in terms of mass values M of 
neighboring nuclei via 3-term formulae, 

2A+ xp (A) = 5 2 M + = 2M(A + 1) - M(A + 2) - M(A) (Al) 

or 

2A cxp (A) = 5 2 M~ = 2M(A - 1) - M(A - 2) - M(A). (A2) 
The 5-term expression is usually considered more accurate, being a half-sum of them, 



A cxp (A) = 5 2 M/2 ee (5 2 M + + 5 2 M~)/4. (A3) 

These simple recipes were used, in particular, in (2- 5] . However, they originate from the 
simplest model of A = const, and the accuracy of such prescription is not a priori obvious. 
To clarify this point we made a calculation which could be considered as a "theoretical 



experiment". We used the GEDF method 28j with the functional DF3 which reproduces 



the mass differences of Eqs. flAlj) . flA2j) type sufficiently well. We calculated first directly 
the right side of Eq. (1A3|) and second, the theoretical gap value with Eq. ( II Op within the 
same GEDF method. The comparison of these two quantities is given in Fig. 7 for the lead 
isotopes and in Fig. 8 for the tin isotopes. We see that for the main part of nuclei under 
consideration the difference between values in two neighboring columns is within 0.1 MeV. 
However, there is several cases where it is of the order (or even exceeds) 0.2 MeV. Leaving 
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FIG. 8: The theoretical predictions for the mass difference ^Mjl (triangles) versus the average 
gap value Ap (squares) for Pb isotopes. Both the quantities are calculated within the GEDF 
method, see the text. 

aside detailed analysis of these "bad" cases we are forced to put a limit of ~ 0.1 — 0.2 MeV 
for the accuracy of the experimental gap determined from Eq. (1A3j) . 
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